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Abstract 

In a one-dimensional periodic nonlinear x^"^^ medium, by choosing a proper mate- 
rial and geometrical parameters of the structure, it is possible to obtain two matching 
conditions for simultaneous generation of second and third harmonics. This leads to 
new diversity of the processes of the resonant three- wave interactions, which are dis- 
cussed within the framework of slowly varying envelope approach. In particular, we 
concentrate on the fractional conversion of the frequency uj {2/3)uj. This phe- 
nomenon occurs by means of intermediate energy transfer to the first harmonic at 
the frequency uj/3 and can be controlled by this mode. By analogy the same medium 
allows "nondirect" second harmonic generation, controlled by the cubic harmonic. 
Propagation of localized pulses in the form two coupled bright solitons on first and 
third harmonics and a dark soliton on the second harmonic is possible. 

PACS numbers: 42.79.Dj 
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1. Introduction 



In recent years a concept of photonic band gap (PBG) materials attracted much attention 
because of their potential applications in many technical and scientific areadlB. The most of 
the theoretical studies of PBG materials focused on the problem of the solution of the dis- 
persion relation that contains an information on the propagation of electromagnetic waves 
through periodic dielectric structures in the linear regime, when the dielectric constant is 
independent of the field intensity. On the other hand, it is well known that the introduction 
of an intensity-dependent refractive index can significantly change the transmission proper- 
ties of a medium. A considerable activity devoted to the investigation of the phenomena 
associated with the propagation of the electromagnetic waves in nonlinear media provided 
a bulk of evidence that PBG structures possessing nonlinearity imply the existence of many 
interesting and useful phenomena which can occur via nonlinear interaction and therefore 
are believed to be equally fruitful as PBG systems based on the linear regime!^. For exam- 
ple, it has been found that the existence of the PBG in which linear electromagnetic effects 
are absent does not exclude the possibility of nonlinear wave propagationiJi. It has been 
demonstrated that for certain values of the input power the radiation with the frequency in 
the stop gap can be transmitted due to the excitation of nonlinear solitary waves which may 
provide practical way of coupling large amounts of optical energy into otherwise unpenetra- 
ble photonic bandgap. If the dielectric has either a positive or negative Kerr coefficient the 
transmission properties of a medium are dramatically changed. For instance, the range of 
the forbidden wavelengths can be altered interactively by the nonlinear medium response 
and this effect can induce an intensity-dependence which has important application in opti- 
cal bistable switching.0 Recently, the study of bistability in periodic x*'^'' materials has been 
reported in Ref. ||. 

On the other hand, the possibility of multiplication or division by an integer of a fre- 
quency of an electromagnetic wave in a nonlinear medium, in particular the second and third 
harmonic generation is of great practical importanceiHi^. In this paper we examine the pos- 
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sibility of the fractional transformation of the frequency which occurs in the photonic crystal 
represented by ID periodic multilayer system possessing x^^^ nonlinearity. This phenomenon 
as well as any scheme of nonlinear conversion must include phase matching mechanism which 
is usually achieved by use of an appropriate birefringent nonlinear crystal or by employing 
some kind of the phase matching which corrects the relative phase between the interacting 
waves by modulating the sign or magnitude of the nonlinear coefficient on the scale of the 
coherence length while the linear part of the refractive index is continuous. In contrast, the 
phase matching in photonic crystals is accomplished by periodic modulation of the linear 
refractive indexffl or by introducing a defect modeEl. The enhanced gain in such structures 
comes about as a result of the interplay of the high electromagnetic mode density density 
in the vicinity of the band edge due to the modified boundary condition£§ and strong con- 
finement of both the pump and higher harmonic signals. Simultaneously, low group velocity 
leads to larger interaction of the interacting waves. The process of the fractional frequency 
conversion which we consider in this paper represents a particular case of three wave interac- 
tion in systems with x*^^^ nonlinearity. Namely, we show that a simultaneous combination of 
u = 3uj ~2uj process and the second harmonic generation in a ID periodic nonlinear medium 
possessing only x*^^) nonlinearity results in the output frequency conversion cu {2/3)uj. 

In the Section II propose ID periodic PBG structure in which the resonant conditions 
for simultaneous SHG and THG and derive a set of coupled-mode equations which allow 
to determine the evolution of the intensities of the interacting modes. In Section III we 
analyze stability of the second harmonic of a constant amplitude in terms of the system 
of the evolution equations and determine the conditions which define regions of stability. 
In Section IV a stationary process is considered in more detail when the initial energy 
is arbitrarily distributed among the modes for phase-matched waves. In Section V we 
identify different types of the evolution of the intensities in terms of the potential governing 
newtonian particle and present the results of numerical simulations which correspond to the 
different regimes. In Section VI a particular solitary wave solution is presented in the case 
when the coefficients 71 and 73 are of the same order. In Section VII we summarize the 



results and discuss both the material the structural parameters which affect the efficiency 
of the processes considered and conditions under which these phenomena can be described 
in the frame of the parabolic approximation. 



2. Statement of the problem 

A. Resonant conditions and phase-matching in ID Photonic Band Gap Structure 

Let us consider a TE polarized plane wave incident on a periodic medium with the dielectric 
perimittivity e{x, lu), e{x, uj) — e{x + L, a;), where L is a period of the structure. We assume 
that the medium possesses a x^^^ nonlinearity which is also periodic with the period L: 
X^'^H^) = X^'^K^ + L). The inclusion of the isotropic second order nonlinearity leads to the 
nonlinear polarization of the form 

PNL{x,t) = x^'\x)E\x,t). (1) 

The Maxwell equations for the TE-polarized wave E = (0, E, 0) propagating along the 
X-axis, i.e. E = E{x, t), lead to the wave equation 

2d^E{x,t) , ^2 r ^^^^ ./Nz:./^ .A J./ . _ 



dx'^ 

where 



+ ^ e{x, t - t')E{x, t') dt' = -An—PNLix, t) (2) 



e(x, t-t') ^ 5(t - t') + ^nxix, t - t') (3) 

and x{x-it) is the linear permittivity [while e{x,uj) is the Fourier transform of e{x,t)\. 

Specifically, we will consider a periodic multilayered structure consisting of two different 
components denoted as a and h and fabricated from the materials possessing x^^^ nonlinearity 
(as the matter of fact it is enough that only one of the slabs is nonlinear). Then we define 
the filling fraction f = a/{a + b) as an important parameter which characterizes 

the stack. 
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Let us assume that the filhng fraction / and the dielectric permitivities of the two types 
of layers e^, e^, are chosen in such a way that there exist three frequencies Uj = uj{qj) 
(j = 1, 2, 3) such that 



where Qj {j = 1,2) are vectors of the reciprocal lattice, are satisfied. In such a system 
one has two resonant conditions for the SHG and third harmonic generation (THG) fulfilled 
simultaneously (if Auj2 = Auj^, = one has the exact phase matching). We notice that one 
of the possibilities to satisfy the mentioned resonant conditions which we call the double 
resonance hereafter, is to choose all uj bordering stop gaps such that (i) all Uj belong to the 
center of the Brillouin zone (BZ) or (ii) uJi and uj^ are at the boundary and 002 is at the center 
of the Bzlil. However, both possibilities are not very interesting from the point of view of 
practical applications, since it is essential to have energy transfer among modes in space, 
while the group velocities of modes bordering a stop gap are zero. Therefore we choose the 
frequencies uj to be shifted towards allowed bands and therefore possessing non-zero group 
velocities Vj. 

We have shown that such resonant conditions can be achieved in the ID periodic structure 
consisting of the alternating slabs of Alo.1Gao.9As and InSb characterized by the dielectric 
constant ecaAs = 10.97 and ejnsb = 16.4 at A = 2/im. Both Alo.1Gao.9As and InSb possess 
the second order nonlinearity with x^'^^ = 1-68 ■ 10^^° m/V and x^"^^ = 1-84 ■ 10~^° m/V, 
respectively. In Fig. |l](a) we present a photonic band structure for ID periodic system 
consisting of the alternating layers of Alo.1Gao.9As and InSb when the filling fraction / = 0.1. 
It is demonstrated that in such a system the resonant conditions given by Eqs. (^) and 
can be satisfied for both wave vectors and frequencies of the fundamental and the final states 
of the second and third harmonics. The band structure shown in Fig. 1(a) is determined 



UJ3 = 3uJi + Acus, UJ2 = luji + Aco'2 



(4) 



where lAcuJ <^ cui and 



gs = 3gi + Qi, q2 = 2qi + Q2, 



(5) 
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by using standard plane wave methods for an ID infinite system and the frequencies are 
expressed in normalized units. In order to correlate our simulations with the experimentally 
accessible systems we choose the lattice constant a = 0.7334/im. The frequencies are taken 
such that the fundamental signal represented by a laser source operating at 2fim corresponds 
to the frequency close to the band gap edge of the third lowest band and the second and 
third harmonic signals with the wavelengths 1/im and 0.667/im, respectively, correspond to 
the frequencies from the vicinity of the 7. and 10. lowest gap, respectively. In Fig. 1(a) we 
indicate the dispersion curve at A = 2/im which corresponds to the fundamental signal by a 
solid line for six lowest bands (as it would be without material dispersion). The broken lines 
refer to the 7. lowest band at A = l/im, which corresponds to the second harmonic and to the 
10. lowest band at A = 0.667/xm which belongs to the third harmonic signal. In evaluating 
the photonic band structure for the fundamental and the final higher harmonic states we 
include the material dispersion in the range of the frequencies considered by substituting the 
tabulated values of the dielectric constantJii at A = 2/im: eAicaAs = 10.97, ejnsb = 16.4; at 
A = 1/im: eAiGaAs = 11-86, ejnSb = 18.23 and at A = O.S/^m: eAiGaAs = 17.77, ej^sb = 13.32. 
The phase matching between the interacting waves was accomplished by varying of the filling 
fraction of the slabs and choosing the components with suitable material parameters that 
constitute a geometrical configuration in which the resonant conditions for both nonlinear 
processes are satisfied. By taking into account the material dispersion of AlGaAs and 
InSb associated with the wavelengths corresponding to the fundamental signal, the second 
and third harmonic we have found that the phase matching between the fundamental and 
the second harmonic and between the fundamental and the third harmonic occurs for the 
opposite propagating waves. In what follows without restriction of generality the geometry 
is taken such that the fundamental (i.e. having the lowest frequency wave) is always forward- 
propagating. In Fig. |l|(b) we show in detail the dispersion curves in the extended scheme in 
the range of the wave vectors and the frequencies in the vicinity of the first lowest band where 
the regions of the curves having positive and negative derivatives correspond to forward- and 
backward-traveling waves, respectively. The phase matching between the forward-traveling 
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fundamental and the second harmonic waves is determined by the intersection of the two 
curves which represent the fundamental signal indicated by the solid line and a dashed curve 
which represents the second harmonic divided by a factor 2. We show that the exact phase 
matching between the fundamental wave and the second harmonic is possible when the 
wave vector q — Qshg — 0.468 and the frequency LoajI'KC — 0.3667, while the exact phase 
matching between the fundamental wave and the third harmonic is possible when the wave 
vector q — Qthg — 0.471 and the frequency ioajl'Kc — 0.3680. The phase matching between 
the forward-traveling fundamental and third harmonic wave is determined by the intersection 
of the two curves which represent the fundamental signal indicated by the solid line and a 
dash-dotted curve which represents the third harmonic divided by a factor 3. Since in the 
system considered the exact phase matching between the fundamental and second harmonic 
and the fundamental and third harmonic occur for the wave vectors qsHG = 0.468 and 
Qthg — 0.471 that are not identical, we introduce parameters Aqshg, ^Qthg, ^^^shg, and 
Acuthg, that characterize the conditions under which approximate phase matching occurs. 
Namely we choose the wave vector q = 0.47 at which the frequency of the fundamental 
signal is uja/2'Kc — 0.3677 and corresponding frequency of the second harmonic divided by 
factor 2 yields LoajI'KC — 0.3667 while the frequency of the third harmonic divided by the 
factor 3 yields loaj^i^c = 0.3680. Then the mismatch in frequency between the fundamental 
and the second the harmonic wave at g = 0.47 is Aa;2 = Aloshg — —0.001 in respect to 
the wave vector q — 0.468 at which an exact phase matching between the fundamental wave 
and the second harmonic occurs. The mismatch in frequency between the fundamental and 
third harmonic wave at q — 0.47 is Aa;3 = Acuthg — 0.0003 in respect to the wave vector 
q — 0.471 at which an exact phase matching between the fundamental wave and the third 
harmonic takes place. 
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B. Evolution equations 



The wave evolution in the processes of higher harmonic generation is well described in terms 
of the so-called envelope function approach!. In the system described above one has to take 
into account three resonant waves, i.e. to express the electric field in the form 

3 

E = J2 Aj4>j{x)e''^^' + c.c. (6) 
i=i 

where (pjix) {j = 1,2,3) are orthogonal and normalized eigenfunction of the eigenvalue 
problem 

[c'^idydx^) + eo(x; cUj)cu|]0j(x) = (7) 

Aj{x) [j = 1, 2, 3) is a slowly varying, compared with (f)j{x)e^'^^* , amplitude of j'-th mode, and 
c.c. stands for the complex conjugate. Using the conventional and well elaborated procedure 
(see e.g.!) one arrives at the system of equations 

^ + + 73^2 Aae^^^--'^-)* + 271^1^26^^-* = (8) 



i + ^^^) + 73AiA3e^(^--^-)* + 7iA?e-^^-* = (9) 



LUs \ dt dx 
where the coefficients are given by 

7i = 27r / -cu, 2cj)0i(x)0i(x)02(a;) rfx (11) 

Jo 

73 = 27r / x^^Ha;; -t^,3t^)0i(x)02(a;)03(a;) rfa; (12) 

JO 

L is a length of the stack, and we assume the following symmetry relations 



20. 



are taken into account and bar denotes the complex conjugation. 

We note that the system of Eqs. (p|)-(p!OD is derived for a finite structure subject to 
periodic boundary condition!. The outcomes, however, can be extended to an infinite 
structure by computing the hmit L ^ oo. The difference between the resuhs obtained 
for finite and infinite structures is estimated to be of 0{L^^). Bearing this in mind in 
what follows analytical results are derived for an infinite (or semi-infinite) structure and for 
numerical estimates of the coefficients 7^ a finite structure are used. Another important 
comment to be made is that we consider also on a semi-line x > 0. Formally to 

do this one should use proper eigenfunctions (t)j{x) of the boundary value problem (rather 
than conventional Bloch functions). For special types of solutions, however, (say possessing 
a symmetry with respect to origin) this leads to a negligible error which very roughly can be 
estimated to be of order of A//, where A is a length of the carrier wave and / is a characteristic 
scale of spatial variations of Aj (in other words to be of order of the small parameter of the 
structure). 

It is worth to point out that in the case of pulse propagation (or more precisely when 
\mit^±oo Aj{x^t) = 0) one of the integrals of the evolution system (|)-(0) is given by 

viWiix) + V2W2{x) + v^W^ix) = (13) 

where 

/oo 
\Aj{x,t)\^ dt. (14) 
-00 

3. Stability of the second harmonic 

One of the solutions of the Eqs. (|])-(|1^) which is of a special importance for the next 
consideration is that of possessing the second harmonic having a constant amplitude the and 
first and third harmonics having zero amplitude. We refer such a solution as background, 
hereafter. Then, considering an infinite system, the question about stability of such a state 
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appears. We consider this problem within the framework of system of the equations (|8|)-([T0D 
for slowly varying amplitudes which are rescaled as follows 

ai = J —Aiexp{-iSuJit), 02 = A2exp{iLpi), as = — A3 exp[i(-(5u;3t+V93+v9i)] 

V W2 i 1 V ""^2 

(15) 

In the case of perfect phase matching {6uj2 = 6003 = 0), we have 
. / dai dai\ _ „ 

? I — + — I + 0203 + 2aia2 = 0, (16) 

/ 1 da2 da2\ _ o / w\ 

.(^-^ + ^2— j+aia3 + a, = 0, (17) 

/ 1 da^ das\ 

I + CX3— + iiaia2 = . (18) 

V W3 Ot J 



Here r = uJiTit, ^ = {uJiTi/vi)x, aj =sign{vj) Wj = \vj\/vi, Tj = \-fj\, ^ = and 
= arg(7j). (Recall that in accordance with our convention vi > 0.) 
We are interested in the stability properties of a solution of the Eqs. ([T6| ) -([18| ) as follows: 
ai,3 = and 02 = pexp[i{uT — kC,)], where u = a2W2k and p is a positive constant playing the 
role of the amplitude of the background. Following the conventional procedure we consider 
perturbations aj of the mentioned solution, i.e. we substitute a„ = a„ exp[r;,i(u;r — + (p)] 
for {n = 1,3) and 02 = (p + 02) exp[2i{ujr — + 0)], where \aj\ <^ p into the Eqs. 
(|T6|)- ([l8|) . Next we linearize the system of equations with respect to the and represent 
Uj = aQj exp[i{QT — K^)] where a^j are constants. Then equations for ai and are singled 
out and we obtain the dispersion relation f{Q,K) = for small excitations against the 
background. The complex roots of this equation indicate instability. After straightforward 
analysis, one obtains the conditions for f{Q, K) = having no complex roots. They read 



Vl 



+ 2 



V3 



< p (20) 
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(fl - V2){V2 - V3,) > 



(21) 



The first consequence comes out from (0): only a background in a form of a plane wave, 
i.e. with /c 7^ 0, is stable. The respective boundary conditions to (p!6|)-(pr8D read 



lim ai^3 = 0, lim 02 = pe 



i(fe^— ui2t) 



(22) 




An important consequence of (|2l|) is that the stable second harmonic can be achieved in a 
structure where phase matching takes place for waves with "ordered" group velocities: either 
Vi > V2 > V3 or vi < V2 < v^. In particular, an unstable case is when the first and third 
harmonics are forward propagating, and the second harmonic propagates in the opposite 
direction. 

As it has been mentioned above in a real configuration it is almost impossible to avoid 
phase mismatches Auj. They, however, result in the same terms of the evolution equations 
as the "rotation" of the background. Hence phase mismatch can either enhance or suppress 
the instability of a static background. Below it will be convenient to consider 002 as a 
reference frequency, and define 6uJi = — Aco'2/2 and 6003 = Au^ — |Au;2 as a mismatch of the 
frequencies of the first and third harmonic, respectively. Then, considering the Eqs. (|TBp- 
(|T8|) subject to the boundary conditions (^), taking into account the phase mismatch, and 
introducing new dependent variables dn = a„exp[i|(A;^ — Wit)], one arrives at the system 
of equations 




(23) 




(24) 




(25) 



where uj 



AiUrj (Tj ~ 

— - — - — jy ■ 




Uj, 61 = ^k{l — W2), and ^3 



3k{w3—W2) 
2w3 
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It follows from the above discussion, that in the particular case when it is possible 
simultaneously satisfy the two conditions 5j = uj the case of quasi-phase matching is math- 
ematically equivalent to the case of exact phase matching with the background in a form of 
a plane wave. Thus, one can associate an effective wave vector of the background with a 
phase mismatch and if this wave vector belongs to the stability region then one can speak 
about stabilization of the background by the phase mismatch. 



4. Stationary process 

Let us consider in more detail a stationary process when A2 does not depend on slow time 
while the temporal dependence of Aj (j = 1,3) is harmonic ~ g-«<5'^j<_ Then the system of 
the evolution equations (|8D-(p!0D can be rewritten in the form 

dai , . 

i-T^ + 0203 + 2aia2 - Viai = (26) 

«cr2— + fliOa + = (27) 

^0-3-^ + 1^0-1^2 - ^3^3 = (28) 
First we concentrate on the case of exact phase matching: Auj2 = Au^ = (z/i = z/3 = 0). 



Then system of the Eqs. (p6|)-(p8|) allows a particular exact solution describing the energy 
transfer from the first and third modes to the second one0. That takes place, however, only 
for a specific relation between the amplitudes of the first and third modes. In the present 
paper we focus on a more generic situation when the input energy is distributed among 
modes arbitrarily. In particular, we are interested in the case when initially the second 
mode is not excited and in the effect of each of modes on the energy transfer to the second 
harmonic. 

To this end we note that the system of the Eqs. (||)-(|l3) possesses two integrals: 

Ar= |ai|' + 2(T2|a2|' + — lasP (29) 

/i 
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[dN/d^ = 0, it substitutes the relation ([I3D in the case of monochromatic waves] and 

H = a\d2 + d\a2 + 010203 + 010203 (30) 

[dH/di = 0). 

For the next consideration it is convenient to separate real and imaginary parts of the 
amplitudes, oi = gi + ipi, 02 = ^2 + io'2P2, 03 = 53 + ia^np^. Then one can consider the 
system of the Eqs. (|6D-(|8D dynamical one generated by the Hamiltonian 



H = (T20-3/igiP2P3 + 2cr2gipiP2 - (^2PiP2q3 + qiq2q3 + (^sl^'P^Ps + qlq2 - p\q2 (31) 

with the canonical Poisson brackets and conjugated variables qj and pj. A remarkable 
property of this system is that the Hamiltonian admits another representation: 

u dq2 . . 

which gives the relation between the amplitude and the phase of the second harmonic. 

In a generic case, when H ^ the dynamics is characterized by permanent energy ex- 
change among modes. In what follows, however, we will be especially interested in situations 
when all energy is concentrated in one or two of higher harmonics. This can happen only if 
H = 0: indeed, if the field of the first or the second harmonic is zero, then H is zero as well. 

The peculiarity of the frequency conversion among modes in periodic media is that 
the modes involved in the process can be either co-propagating or counter-propagating. 
Naturally, the behavior of the system strongly depends on the case under consideration, as 
we already have seen on example of the stability. We can distinguish four different cases: 
(i) forward propagating waves {a2 = 0-3 = 1), (ii) backward propagating second harmonic 
{(72 = — 1 and o"3 = 1), (iii) backward propagating third harmonic {a2 = 1 and = —1), 
and (iv) backward propagating second and third harmonics (cr2 = cr3 = —1). Notice that 
the last case occurs in the structure depicted in Fig. 0. 

In order to study the dynamical system for qj and pj we observe that it follows from ( |52[ ) 
that when H = one obtains 52/^2 = cot02 where 62 is the phase of the second harmonic 
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(hereafter we speak about the phases of Oj) which does not depend on the ^ and thus is an 
integral of motion. On the other hand, by computing p2 from (BTI) we obtain 



qiQs + fj'PiPs + qi-Pi 

The last formula allows one to express the functions qi, pi, q^, and ps through the other 
ones. 

The next consideration is devoted to the case when the second harmonic has phase shift 
7t/2 with respect to the first and third ones: 

02 = ±n/2. (34) 

This means that q2 = and p2 7^ at any point in space. Then by introducing a new 
variable 

Xp = a2 I^V2{0di' (35) 

one can linearize the dynamical system which now takes the form 
dqi 



dXp 



^ - gs + 2gi = (36) 



Z| + = (37) 

# - ^ = (38) 

dXp qi 

The first two equations are linear and thus are trivially resolved with respect to qi and 

= Cie''+^^ + Cse'"-^^' (39) 

gg = _Cicj_e'^+^^ - Caw+e'"-^^ (40) 
where 
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1 ± \/l — cr-^jj, if a-iji < 1 



(41) 



— 1 ± iy/a^fi — 1 if asjj, > 1 



Then for pi we obtain 



Pi = Poie' 




(42) 



The constants Ci, C2, and poi are determined by the amphtude of the field at the input of 
the structure. 

5. Fractional frequency conversion 



at ^ ^ 00. In the case when they exist, stationary points correspond to situation where 
the energy of the electromagnetic field is either distributed among modes or concentrated 
in one of them and energy exchange among modes does not occurs any more. Then one 
can distinguish two cases: (i) Xp ±00 what happens when ±(72^2**^ > 0, and (ii) Xp — >■ 
X^'^*) = const, what happens when = 0, where P2^*^ = lim^_,oo (0 • The first situation 
corresponds to the case when the energy is converted to the second harmonic and the second 
one corresponds to the case when the energy output of the second harmonic is zero. In this 
section we concentrate on the former case. 

It follows from the Eqs. (|39D, (|40|), and (^) that in the case of forward propagating third 
harmonic (cts = 1) a bounded nontrivial solution can exists only if poi = 0. Moreover, a 
convergent solution of the Eqs. (0), (|37|) can exist only when a2P2^^ > what corresponds 
to the Qi'^^ = q^^^^ =0. In physical terms this conclusion means that the energy initially 
distributed among three modes is totally transferred during the process into the second 
harmonic. The necessary condition for such process requires the first, ai, and third 03, 
harmonics to be forward propagating waves and have phase difference either or tt while 
the second one, 02 has a relative phase shift ±7r/2 [the actual phases of the modes are defined 
by the last conditions and the relations ([T5|)]. 



Let us look for stationary points of the Eqs. (^6])-(^), i.e. for solutions tending to constant 



15 



In order to understand which parameters of the input signal result in stationary points 
we observe that Xp as a function of ^ can be treated as a coordinate of an effective newtonian 
particle governed by the equation (the case as = —1 is included in the next consideration, 
as well) 

The potential U{Xp) is easily computed from the equation for q2 and (^), (|40| ) in terms of 
elementary functions for arbitrary value of /i. It has, however, a rather cumbersome form 
to be represented here. Instead, we mention that the dynamics of the effective particle 
essentially depends on whether (a) < a^n < 1, (b) a^fi > 1 or (c) a^fi < and /i > 1. 
In order to simplify the discussion of these cases we assume that there is no input second 
harmonic. This means that initially (i.e. at ^ = 0) Xp = 0, dXp/d^ = i.e. the energy 
of the effective particle is determined by the point of intersection of the energy curve with 
[/-axis, and that the type of the motion depends on the value of the potential energy at 
Xp = 0. In Fig. we display three typical configurations of U (Xp) which correspond to the 
above cases. They represent different kinds of motion: Fig. ^(a) corresponds to motion of 
the particle to the right which ends up in the steady state motion. This type of motion 

(st) 

means that Xp — s ^ind thus corresponds to the total energy transfer from the first 
and third harmonics to the second one. In Fig. Q(b) the particle is initially in a potential well 
and thus undergoes oscillations around the local minimum of the effective potential. This 
type of motion, corresponds to periodic energy exchange among modes. Finally, Fig. H(c) 
corresponds to exponentially growing solution which indicates the instability within the 
framework of the parabolic approximation. 

The first two types of the behavior described above are represented in Figs, ^(a) and|^(b), 
where we show the results of the numerical simulations obtained by solving the dynamical 
equations (^)-(^) for slowly varying amplitudes Aj (j = 1, 2, 3). We display the evolution of 
the intensities \Aj\'^ {j = 1, 2, 3) as a function on ^. In Figs. 0(a), (b) we have chosen the same 
initial values of the amplitudes of the input signal and the same structure parameters as those 
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used in Figs, ^(a) and @(b). In Fig. ^(a) energy transfer from the first and the third harmonic 
to the second one is shown. The time of the energy transfer is inversely proportional to the 
square root of the field intensity which is a direct consequence of the evolution equations. 
Taking into account that the amplitude of the first harmonic is relatively small, this situation 
can be viewed as frequency down-conversion - 003 ^ {2/ 3)003 - in the case when the input 
energy is dominated by the pump signal with the frequency u!{= u^). Then after the 
transition time the energy is transferred into the mode with the fractional frequency |u; 
(= UJ2). The fractional down-conversion is observed even when the first harmonic is negligibly 
small (but nonzero) at the input. The existence of the first harmonic in this process is 
fundamental. Indeed, = const., Ai = 0, and ^2 = is a solution of the dynamical 
system (§)-(|10|) and thus the first harmonic is necessary in order to initiate the process of 



the down conversion. Moreover, at small enough intensities of the first harmonic we have 
observed local amplification of the first harmonic. 

The case which corresponds to the periodic exchange of the energy among modes is 
shown in Fig. ^(b). Here we are close to the minimum of the potential energy and the 
interaction among modes is rather weak (the second harmonic is not zero but has a very 
small amplitude). Almost the whole energy of the electric field is concentrated in the third 
harmonic. However it is possible to find parameters at which interactions among modes are 
strong and even locally the whole energy is concentrated in one or two modes. 

As it is evident, by a proper choice of initial parameters in the case of a "potential " 
depicted in Fig. ^(b) an effective particle can be placed outside the potential well. Then 
the fractional frequency conversion occurs. This situation is given in Fig. 0(c). Finally, in 
Fig. |](d) it is shown the exchange of the energy among modes in the presence of mismatch. 
Although the second harmonic is still generated (i.e. fractional down conversion occurs) at 
specific distances which are placed almost periodically one observes enhancement of the first 
harmonic. 

Returning to the example of a real structure described in Section 2 one cannot provide 
frequency conversion to one of the frequencies because of strong instability (this instability 
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however is an artifact of the parabohc approximation: direct numerical simulations of this 
case will be reported elsewhere). Meantime, when the third harmonic has na amplitude 
much higher than the first and second harmonic, one can achieve a stable propagation of 
the modes, which, however, manifest very weak interaction. This situation is shown in 
Fig. The stability of the picture represented is very sensitive to the relative phases and 
the amplitudes of the waves. 



6. One solitary wave solution 

Let us now return to evolution equations (p3D -(p5D. As it has been mentioned, the phe- 
nomenon described above is essentially based on the geometry of the structure. As a conse- 
quence the effective nonlinear coefficients 71, and 73 can be of different orders. In particular, 
one can achieve the relation |7i| ^ I73I. Then the equations (^)-(^) in the absence of phase 
mismatch are reduced to the conventional system for the SHG which is integrable by means 
of the inverse scattering technique^. Another integrable case occurs when I73I ^ |7i|. Then 



the Eqs. (|23D-(|25D are reduced to the system describing decay of the pump wave, which is 
the third harmonic, and the back process: A3 ^ Ai + A2. A characteristic feature of the 
last process is that the frequency of the three waves are related by @). 

In the case when 71 and 73 are of the same order it is still possible to find a particular 
solitary wave solution having the form of the coupled bright solitons of the first and third 
harmonics and a dark soliton associated with the second harmonic. It reads 



^a±ui \ (3ui ' ^' ^' j 



52 = "^^7^ I + tanh(/5C) 1 (45) 



^ia\pu\ 1 

"3 = 1 2 W 

1 — jjLa%.ui cosh(p(, j 

Here 



- Vt)v2 



V2 - VVi 



V is a velocity of the pulse, 



u — 



V3(T3 - VVi 
(1 - V)V3 ' 



VVi 

ui = as 

^^3 




fj,{fj, + u) ± \/u{u — /i) 



liui{ii + 3u) 



and the phases of the frequencies of the modes must be chosen to satisfy the relation 

^ ^ (1 - fialuiY 
Di a^uidia^ui — 2) 

7. Discussion and Conclusion 

Within the slowly varying amplitude approximation we derived a set of coupled mode equa- 
tions which describe the evolution of the intensities of the resonant waves. We used these 
equations to study double resonant processes in x^^^ nonlinear periodic media in which two 
matching conditions for the generation of the second and third harmonic are simultaneously 
satisfied. We show that such conditions can be achieved in one-dimensional photonic band 
gap structure by choosing a proper combination of geometrical and material parameters. 

Speaking about possibilities of experimental fabrication of structures with double phase 
matching we notice that it is more difficult to match parameters, than say for the second 
harmonic generation only. In particular it is intriguing when actual experimental values of 
the refractive index are taken into account. Meantime it seems that there are no restriction 
for possibility to find a proper geometry almost for any x^^^ material. Another possible 
technological difficulty could be a lattice mismatch of the layers. In this context it is worth 
noting that a periodic structure is rather fiexible in the sense that it allows inclusion of 
new components which being much thinner than the main nonlinear slabs do not affect the 
phenomenon of frequency conversion. Namely in the case when chosen layers possess large 
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lattice mismatch, slabs of the third kind having width much smaller than a or 6 could be 
included between nonlinear layers. 

It has been shown that the use of double resonances allows one to obtain difference 
frequency generation. A particular example is ^ (l)"^- This however is not direct 
conversion but it occurs with participation of the mode with the frequency cj/S. 

Parabolic approximation which for a given scaling results in a dynamical system is gen- 
erated by Hamiltonian H. The non-zero Hamiltonian corresponds to the permanent energy 
exchange among the modes, while H = corresponds to the case when the energy can be 
(subject to definite conditions) concentrated in one or two higher harmonics (that is the 
case of fractional conversion). In processes of the frequency conversion in periodic media 
the modes can be either co-propagating or counter-propagating and consequently the be- 
havior observed is very different. In particular, one can distinguish the following types of 
the dynamics: (a) the total energy is transferred from the first and third harmonic into 
the second one, (b) a periodic exchange between the modes occurs and (c) exponentially 
growing solution, which represents instability within the parabolic approximation. The qual- 
itative picture based on the analysis of the motion of the effective particle in the potential 
which indicates different types regimes is very well reflected by the results obtained from 
the numerical simulation carried out to solve dynamical equations. 

The evolution equations governing three-wave processes in the presence of double reso- 
nances allow a solitary wave solution in a form of coupled two bright solitons (they corre- 
spond to the first and third harmonics) and a dark soliton on the second harmonic. 

Finally, we also notice that periodic structures, involving geometrical factors rather than 
only material properties can provide us with a large diversity of matching conditions. In 
particular one can obtain double resonances as follows [c.f. the Eqs. (^) and (^] 

iU3 = 2tU2 + Acua, iU2 = 2uji + Alu2 , (48) 

and 

gs = 2g2 + Qi, q2 = 2gi + Q2 , (49) 
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which result is a different kind of three wave interactions. Sohtonic solutions generated by 
the respective interactions in the presence of dispersion have been recently considered i ni. 
Naturally, the dispersion can lead also to solitons in the case of double resonances defined 
in the Eqs. (§, (§. 
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FIGURES 

Fig. 1. (a) Photonic band structure of ID periodic structure consisting of alternating slabs 
of Alo.1Gao.9As with Ca = 10.97 and InSb with = 16.4 at A = 2//m. (solid curve) in which 
fractional frequency conversion takes place. The broken curves refer to the 7. lowest band at 
A = Ifim which corresponds to the second harmonic and the 10. lowest band at A = 0.667/Ltm 
which corresponds to the third harmonic signal, q is measured in the units 27r/(a + 6)). (b) The 
detailed picture of the region of the photonic band structure shown in Fig. la in extended zone 
scheme in which both frequencies and the wave vectors that satisfy the resonant conditions for 
simultaneous SHG and THG. The solid curve indicates the calculated dispersion curve near 2nm, 
the dashed line refers to the region near l/xm with values divided by a factor 2 and the dash-dotted 
line refers to the region near 0.667/xm with values divided by a factor 3. The exact phase matching 
between the forward-travelling fundamental and oppositely travelling second harmonic occurs when 
Q = QSHG = 0.468, while the exact phase matching between the fundamental wave and third 
harmonic is possible for oppositely propagating waves when q = qrHG = 0.471. q is measured in 
the units 27r/(o + 6)). 

Fig. 2. Examples of the effective potential U{Xp) for the structure parameters and initial ampli- 
tudes of waves as follows: 71 = 82.14 + 78.48z, 73 = -14.57 - 5.07i, (a) vi = 0.5382, V2 = 0.2808, 
V3 = 0.2433, Ai = 0.13 • lO'' V/m, A3 = 0.31 • 10^ V/m, A2 = 0, fi ^ 0.1225; (b) vi = 2.382, 
V2 = 0.2808, vs = 0.06433, Ai = 0.1 • 10*^ V/m, A3 = -0.9 • lO'^ V/m, A2 = 0, ^ ^ 2.051; (c) 
vi = 0.5382, V2 = 0.2808, V3 = -0.2433, Ai = 0.13 • 10^ V/m, A3 = 0.31 • 10^ V/m, A2 = 0, 
H ^ 0.1225. 
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Fig. 3. Evolution of intensities of the first (solid line), the second- (broken line) and 
third-harmonic (dashed line) signals, (a) the same structure parameters and input amplitudes 
as in Figs. 2a; (b) the same structure parameters and input amplitudes as in Figs. 2b; (c) the 
same structure parameters as in Fig. 2b and the input amplitudes Ai = 0.1-10^ V/m, ^3 = 0.4- 10^ 
V/m, A2 = 0; (d) the same structure parameters and input amplitudes as in Fig. 2a except the 
phase mismatch Auj2 ~ —0.002 and Auj^, ~ 0.0009 was included. The intensities represented are 
normalized to [x^^^]^^ and time is measured in uj^^ units. 

Fig. 4. Evolution of intensities of the first (solid line), the second- (broken line) and 
third-harmonic (dashed line) signals in a structure depicted in Fig. where the input amplitudes 
Ai = 0.255 • 10^ V/m, ^3 = -1.996 • 10^ V/m, A2 = 0. 
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